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Abstract 

An analysis on the damped quantum search by exploring the rate at 
which the target state is obtained. The results were compared with that 
of the classical search since the standard Grover's algorithm does not give 
a convergent result if the number of target state is unknown. For a large 
number of target states, the classical and the damped quantum search 
give a similar result. However, for intermediate values of the target size 
the damped quantum search gives a higher probability of success than the 
classical search. Furthermore, we also made an analysis on the average 
number of iterations needed to obtain at least one of the target states. As 
the number of target states is reduced, the damped quantum search gives 
a better result than the classical search. The results coincide if the size of 
target state is comparable to the size of the sample. 

1 Introduction 

Harnessing the power of quantum mechanics in computation promises a new 
era in computer science. Two well known quantum algorithms that shows su- 
periority over their classical counterparts are due to Shor and Grover. Shor 
demonstrated that the problem of finding the prime factors of an integer can be 
solved efficiently on a quantum computer pQ. Grover showed that the problem 
of conducting a search through some unsorted database can be sped up on a 
quantum computer [2]. While Grover's algorithm did not provide an astonishing 
speed up as Shor's algorithm, the widespread applicability of the search-based 
problem has excited significant interest in Grover's algorithm [3J. 

We can describe the search problem as follows: Consider a database with N 
unsorted objects, among which M of them are the desired objects, how many 
times on the average do we need to search in order to find the first desired 
object? Using classical computation it will take N/M search on the average 
before the first desired object is obtained. A detailed derivation of this result is 
in [3] . If we use the quantum search algorithm it will take */ N/M on the average 
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before we obtained the first desired object [3] . Experimental implementations of 
Grover's algorithm on a small scale quantum computer are reported in [6] and 
they were able to show that the search requires fewer steps than on a classical 
computer. 

The introduction of the Grover's search algorithm spark the interest of some 
researchers to address its shortcomings [U [7J [5J E]. It has been established 
already that the Grover's algorithm provides a quadratic speed up over its clas- 
sical counterpart and was shown to be the best oracle based search algorithm 
[TU] . However, the Grover's algorithm is not without limitations. One of these is 
knowing the number of target states before the procedure is performed. With- 
out such knowledge, the process simply oscillates. One way to approach this 
dilemma is to introduce a damping parameter to suppress the mentioned os- 
cillation. This was done by Mizel in [11] where he attached an external spin 
with the corresponding damping parameter. He demonstrated that there exists 
a critical damping which separates the classical result from quadratic speed up 
of Grover's result. 

The quantum search algorithm [5] starts from an initial state vector to a 
target state with a number of iterations. This leads to a quadratic speed up 
over its corresponding classical algorithm in searching an unsorted database. 
To have an optimum result, one needs to know the number of target states. 
Boyer et al. provided a tight analysis of Grover's quantum search algorithm 
by giving a simple closed-form formula for the probability of success after any 
given number of iterations |10| : 



where k and I are the probability amplitudes of the target and nontarget states 
respectively. Here M is the number of the target states, N is the total number 
of items in our database, j is the number of iterations, and the angle 6 is defined 
so that sin 2 6 — M/N. This allows us to determine the number of iterations 
necessary to achieve with almost certainty of finding the target state, as well as 
the upper bound on the probability of failure. The closed-form formula for fc» 
and lj could leisurely be derived by mathematical induction. 

Mizel proposed a quatum search algorithm that utilizes dissipation in order 
to make the search robust even without the knowledge of the number of target 
states [11] . However, providing a closed- form formula for the case of damped 
quantum search is not trivial. 

In this work, we apply the damped quantum search in an unsorted database 
of 8- and 12-spin Ising systems and provide an analysis of the damped quantum 
search. 



k 



= ^-sin((2i + l)fl) 
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2 Damped Quantum Search in an Ising Spin 
System 

The recent quantum search algorithm introduced by Mizel is appropriate in 
searching an unsorted database without the prior knowledge of the number of 
target items M. We expound the significance of Mizel's algorithm by applying 
it to a one-dimensional Ising spin system database. We search for an eigenstate 
with 8 and 12 spins using the damped quantum search algorithm. 

An n-spin Ising system models a ferromagnet with n magnetic dipole mo- 
ments. Its Hamiltonian is given by H = —eY^i=i s i s i+ii where e is the interac- 
tion energy. It can be shown that the Hamiltonian is a 2 n x 2™ diagonal matrix. 
Here, s is the cr z -Pauli spin matrix defined in a two-dimensional Hilbert space 
spanned by | f) and | \), representing the spin- up and spin-down states. Thus, 
we can construct an oracle that behaves as follows: 

[-1 ifHu = X 
[+1 ifH i}i J=\ 

where A is the desired energy eigenvalue. 

We seek on 8- and 12-spin Ising systems for at least one of their eigenstates 
of a certain eigenvalue by applying the damped quantum search. For compari- 
son, we also show the classical result and the quantum search result assuming 
ignorance of the degeneracy M. There are N = 2 8 = 256 distinct configurations 
that we take as items in the database of the 8 spins system, and N = 4096 for the 
12 spins. Our task is to locate at least one of the eigenstates that corresponds 
to a certain eigenvalue, A, for n = 8 and n = 12. 



2.1 Numerical Computation 



To implement the damped quantum search algorithm in an n-spin Ising system, 
we find the degeneracy M of the eigenvalue corresponding to A for validation 
purposes. It would be convenient for computational intent to write the Hamil- 
tonian as follows: 



H = 



71 — 3 



S ® /®(™- 2 ) + Y] J®« j®(n-t-2) + j®(n-2) ^ g 



where S = s <g) s and / = J® 1 is a 2 x 2 identity matrix, 





S 



-1 



I 



1 
1 



(4) 



(5) 



Since s and I are all diagonal matrices, then H is also diagonal and finding the 
eigenvalues is straightforward. 

Let us now implement the quantum search algorithm. We construct the 
initial state and the Grover rotation operator G using our knowledge of 
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N = 2" and the degeneracy of the desired A, M 
N - 2M 



G = 



N 



NM - M 2 



N 2 



-1 



(6) 



The probability that the target state has been reached were recorded in each 
iteration. 

On the other hand, for the damped quantum search we assume ignorance of 
the M. We start with the construction of the square matrix 



Tr(p'X) 
Tr(p'Z) 
Tr(p') 



cos 29 cos 4> 
— sin 26 cos < 
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where p = \ip){ip\, X = 

eter cf> is defined as cos < 
the column vector 
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1 



and Z 





-1 



Tr(pX) 
Tr(pZ) 
Tr(p) 

(7) 

The damping param- 



= (1 — sin6')/(l + sin 9). Equation[7] will be applied on 

(8) 



Tr(pX) 
Tr(pZ) | = 
Tr(p) 

The value of the third entry of the resulting column vector is then subtracted 
from one and the result is plotted as a function of iteration. 




2.2 Probability of Success 

Figure Q] shows the behavior of the probability of finding the target state for a 
given M in an 8-spin Ising system using the Grover's quantum search, classical 
search and damped quantum search as the number of queries increases. 

The oscillations of the probability is evident in Grover's quantum search. If 
we do not know the degeneracy of the desired eigenvalue, then the search is not 
optimized. Furthermore, in the damped quantum search, the probability that 
the target state will be found and the external spin flipped steadily approach 
unity as we increase the number of queries. This is similar to the classical 
limiting case. The advantage of the damped quantum search is that it preserves 
the 

signature 0(y / M/N) number of queries [11]. Consider for example Figure Q] 
(a), for n = 12 with M — 22, there are 10 iterations needed to make Grover's 
search robust. This gives a probability of 99.9%. The damped quantum search 
locates one of the target states after 32 queries with a probability of 99.5% as 
shown in Figure [T] (b). Ignorance of the degeneracy requires additional queries 
but are minimal compared to the classical search. The latter takes 247 queries 
with the same probability. 
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Figure 1: (a) and (c): The oscillations of the probability of finding the target 
state using quantum search, (b) and (d): the increase in the probability that 
the final system will collapse to ® | t)(t I using damped quantum search 

and classical search for N — 4096 at a given M. 
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Table 1: Minimum expected number of iterations before success for classical 
search and damped quantum search. We take N = 256 for 8 spins and assume 
ignorance of the degeneracy M for a given energy A. 



A 


M 




^dqsmin 


-^csmin / -^dqsmin 


±7e 


2 


39.1796 


19.1233 


2.04879 


±5e 


14 


7.4034 


6.5365 


1.13260 


±3e 


42 


3.2121 


3.2833 


0.97831 


±le 


70 


2.3507 


2.3986 


0.98003 



Table 2: Minimum expected number of iterations before success for classical 
search and damped quantum search. We take N = 4096 for 12 spins and 
assume ignorance of the degeneracy M for a given energy A. 



A 


M 


Ecsmin 


Edqsiniri 


^csmin / -^dqsmin 


±lle 


2 


539.9602 


79.2205 


6.8159 


±9e 


22 


55.1413 


23.2660 


2.3700 


±7e 


110 


13.2779 


9.8111 


1.3533 


±5e 


330 


5.5076 


5.1825 


1.0627 


±3e 


660 


3.2537 


3.3259 


0.9784 


±le 


924 


2.6085 


2.6638 


0.9792 



3 On the Average Number of Iterations 

To generalize our results, we employ the concept of the average number of 
iterations. We start with j iterations. After the j th iteration, we check if 
the target state has been found. If the result is positive, we stop the search. 
Otherwise, in case of failure we restart the iterations. After the first j iterations 
has failed, the average number of iterations, E(j), is given by 

where P{j) is the probability that the target has been found after the j th iter- 
ation. In line with this, we want to obtain a certain value of j that will give 
the minimum average number of iterations, E csm i n for the classical search and 
Edqsmin for the damped quantum search, before success. If a closed form of 
P{j) is available, then we calculate E'(j) and equate it to zero. This is not the 
case for the damped quantum search since the closed form of P(J) cannot be 
obtained easily and the form of P(j) varies for each iteration. Therefore, we 
can plot E(J) with respect to j using [H] to obtain the minimum graphically . A 
summary of all minimum average number of iterations before success for 8- and 
12-spin Ising system are summarized in Tables [T] and [2j respectively. The over- 
head in the number of iterations of the classical search relative to the damped 
quantum search is obtained from their ratio. 

Using the energy eigenvalues as the search marker of the target states makes 
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Figure 2: Expected number of iterations before success of classical search (+) 
and damped quantum search (o) as a function of the number of iterations for 
N = 256 at a given M. 
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Figure 3: Expected number of iterations before success of classical search (+) 
and damped quantum search (o) as a function of the number of iterations for 
N = 4096 at a given M. 
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the search process dependent on the corresponding degeneracies. The behav- 
ior of the average number of iterations is the same for both the classical and 
the quantum search. However, the difference of the classical and the quantum 
search is evident when the target state is one of the highly excited state where 
the degeneracy is very small. The ratio of the average number of iterations 
approaches unity as the energy eigenvalues, which serve as the marker, ap- 
proaches the ground state energy. This shows that the quantum search process 
approaches the classical result if there are several states with the same marker. 

4 Conclusion 

We have shown the advantage of damping the quantum search. Application of 
the Grover's quantum search in an Ising system reveals its oscillatory nature that 
leads to a frail search. This dilemma becomes more evident if we are unaware 
of the degeneracy M in advance. The damped quantum search provides a 
fixed point of convergence, that is, increasing the number of queries will get us 
closer to the target state. Furthermore, for large number of target states, the 
damped quantum search transitioned to classical search. We consider this as an 
advantage in quantum searching because the number of iterations needed by the 
search algorithm increases with M for M > N/2, thus we prefer using classical 
search in this case. On the other hand, the overhead on the minimum number 
of iterations before success of the damped quantum search over the classical 
search is apparent for small number of target states. 
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